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Abstract 

The general Helfrich shape equation determined by minimizing the cur- 
vature free energy describes the equilibrium shapes of the axisymmetric 
lipid bilayer vesicles in different conditions. It is a non-linear differen- 
tial equation with variable coefficients. In this letter, by analyzing the 
unique property of the solution, we change this shape equation into 
a system of the two differential equations. One of them is a linear 
differential equation. This equation system contains all of the known 
rigorous solutions of the general shape equation. And the more general 
constraint conditions are found for the solution of the general shape 
equation. 
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In recent years much attention has been payed on the investigations of the equi- 
librium shape of lipid bilayer vesicles in aqueous solution by physicists theoretically 
and experimentally [1-10]. The amphiphilic molecules such as phospholipids assem- 
ble to form bilayers in water and at low concentration close to take single shells 



which are called vesicles [11]. These bilayer vesicles are considered as the models 



for biomembranes or cells ||12|| . The elastic property of the phospholipid vesicles is 
related to their behaviors about their formation, stability, size, fusion and budding 
processes. The Helfrich theory [|l| and general shape equation || play an important 
role in the study of the equilibrium shapes of lipid bilayer vesicles. The equilibrium 
shape is determined by minimizing the curvature free energy [[TJ 

F = ^-<f(c 1 + c 2 -c ) 2 dA + \JdA + ApjdV (1) 

where k c , ci, c 2 and c are bending rigidity, two principal curvatures and spontaneous 
curvature, respectively. A and Ap are Lagrange multipliers due to the constraints of 
constant volume and area and the former denote a tensile stress and the later is the 
osmotic pressure difference between outer and inner media. The general equation 
of mechanical equilibrium for vesicle membranes is a non-linear differential equation 
with variable coefficients, which is derived from the above curvature free energy. As 



the view in Ref. fl3|| , the general shape equation of the axisymmetric vesicles is a 
nonlinear third-order differential equation and has the form: 
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where ip is the angle between the surface normal and the symmetric axis, and p is 
the distance from the symmetric axis. 



The very interesting work by Naito, Okuda and Ou-Yang |T3| presented a new 
exact solution of general shape equation of axisymmetric vesicles, which enriches 
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the solution family of the shape equation. Now we have known that this general 
shape equation has the solutions of spheres ||, circular cylinders, Clifford torus f|, 
Delaunay's surfaces, nodoidlike shapes and unduloidlike shapes [fL3"l| . These solutions 
were obtained by guessing and their correctness were checked by taking them into the 
shape equation. Are there any method to solve this nonlinear third order differential 
equation for finding its solutions? It is just one of the aims in this letter. And we 
find that all the above shapes of vesicles are the particular solutions of a linear 
second order differential equation. The more general constraint conditions exist for 



the "three-term" solution comparing with the one in Ref . |L3| . In the following, we 
discuss them in details. 



Setting 

\f = sin^(p), 
we can write down the general shape equation as 
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The power series solutions for vesicles were studied in the interesting work |TJJ where 
they expanded the solution of the shape equation as the infinite sum of a power 
series for the distance from the rotational z-axis. In order to discuss the property of 
the cutoff of the power series solution we now put the solution of the general shape 
equation into the form 

M 

* = E *sp s - (5) 



By an analysis of the shape equation we decide the parameters m and M in the 
above expression as follows. The terms related coefficient a 3 (\s \ > 1) are 
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in the shape equation (4). It gives that 

0-l)(s + l)(3s-l) = 0. (7) 

This means that M — 1, m — — 1 and expression (5) is 

^ = — + a + a x p. (8) 
P 

We call it as the "three-term" solution of the general shape equation. On the side, 
from the Ref. J7|, we know that the shape equation (4) has the first integral form. 
So it can be written as 

W ~ l )p^~ \ R p\^f + ^ 2r2 p 2 - - c oR 2 p+(^ + §)Rp = f(ph (9) 

where R = \& / p and the function /(p) satisfies the differential equation 

>» + >> + £ = °. do) 

which gives that 

2 k p 

where C is a integral constant. In this way, the general shape equation becomes as 
a nonlinear second-order differential equation: 

,,,, 9 . gPF Y /dY\2 dY /A Cn\,, 9 An , „ 9 , . 

< y - "V - 2 ( dp ) + - c » y + (* + i) y " + -i" ~ Cp = (12) 

with the use of the notation Y = p 2 R = p^. By considering the discussions in (6-8), 
that is, M — 1, m — — 1, and the known solutions given in Refs. || |7|, IT3, we make 
the Ansatz: 

where Ci, C*2 and C3 are the arbitrary constants. Then, from equation (12), we have 



that 
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It is a linear differential equation of the second order and can be changed into the 
first order one: 

f ~-Y = ^ + ^ + C 3 , (15) 
dp p p 2 p 

dY 

where Y = ——. Then we get solution 

dp 

Y = \{C 3 p 2 - d) Inp + l(Ax - ^)p 2 - C 2 p + A 2 (16) 

where A 1: A 2 are the integration constants. By substituting 

dY „ , „ ^ Ci 
— = G 3 p In p + A x p - C 2 - — , 

ap 2p 

_ = C 3,„ P + C 3 + ,4 1 + ^, 

and relation (16) into equation (13), we can get the constraint conditions of the 
physical parameters A, Ap, c and the constants Ci,C 2 ,C 3 ,A 1 ,A 2 : 

2p 4 (lnp) 2 C 3 2 (2c - C 3 ) + p 4 lnpC 3 (8A lCo - 4^ - 4c 2 - 4c C 3 + 3C 3 2 - y) 

8AA , „ ^ 2 ^ 3 , 4C 3 A 
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|+4p 3 lnpC 2 C 3 (-4co+3C 3 )+4p 3 C , 2 (-4A 1 co+3A 1 C 3 +2c 2 +2c C , 3 -2C : 
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-4Ai A 2 C 3 + 2A!C 2 2 - 4A 2 c 2 - 4A 2 c C 3 + 4A 2 C 3 2 - + 8C + 8c C 2 - 9C 2 2 C 3 + 8C 3 ) 



+ 16p lnpA 2 C 2 C 3 + %pC 2 {2A x A 2 - 4A 2 c + AA 2 C 3 - C 2 + 2) - 16 lnpA 2 C 3 
+ %A 2 {-2A X A 2 + 2A 2 c - 2A 2 C 3 + C 2 2 ) = 0. (d = 0) (17) 

In this way, expression (16) gives the logarithm solution 

Y = c p 2 ln(p/p ), (18) 

when we choose that 

X = Ap = C 1 = C 2 = A 2 = 0, C 3 = 2c , 

C = -2c , Ai = -2c In po + c . (19) 



When d = C 3 = and 

C = \a\A 2 - A 1 A 2 c - -A X C\ + A 2 (| + - c Cl (20) 



2 1 ' " 4 ^2 A; 

we obtained the solution 



r = ^p 2 -C 2 p + A 2 (21) 



under the conditions of 



9 9 2AxA 2Ap . , 

A\c - A x c\ - — ^ ^ = 0, (22) 

^ 2 (^ + y-2A lC o) =0, (23) 

C , 2 (2v4 1 A 2 - 4A 2 c - C 2 2 + 2) = 0, (24) 

A 2 (2A l A 2 - 2A 2 c - C\) = 0. (25) 

This solution contains the one presented in Ref. [0 due to the above conditions 
have the more general form. It degenerates to sphere, circular cylinders, tori and 
Delaunay's surface in the limiting cases. For instance, in the case of A 2 = C 2 = 0, 
it reduces to the sphere solution with the condition AprQ + 2Ar — kc (2 — c r ) = 
where r = 2/A\. And it gives the Delaunay's surface for C 2 = and A 2 ^ under 
the condition cq\ + Ap = 0. 

As the conclusions, on the basis of the discussion of the unique property for the 
solution with finite terms about the general Helfrich shape equation, we change the 
general shape equation into a linear and a nonlinear differential equations. This 
equation system contains all known rigorous solutions of the shape equation. We 
find also that the restricted conditions (22-25) of the polynomial solution (21). These 
conditions have the more general form comparing the known solutions ||13|| . These 
known solutions are all the particular solutions of a linear differential equation of 
the second orders. And the above discussions remind us that, except the solutions 
presented in this letter, the other exact solutions which we guess it exists for the 
shape equation can not be denoted by the elementary functions. It is a challenged 
problem to find out them with the proper choice of the Ansatz. Another interesting 
problem is about the stability considerations of the new solution as in Ref. [IB]. We 
hope the above discussions will be helpful for the study of biomembranes. 
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